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ABSTRACT 

We present the results of a search for extraterrestrial electron antineutrinos (i^e's) 
in the energy range 8.3 MeV < 3^^ < 31.8 MeV using the KamLAND detector. In an 
exposure of 4.53 kton-year, we identify 25 candidate events. All of the candidate events 
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can be attributed to background, most importantly neutral current atmospheric neu- 
trino interactions, setting an upper limit on the probability of solar z^g's converting 
into z^e's at 5.3 x 10~^ (90% CL), if we assume an undistorted Vf, shape. This limit 
corresponds to a solar V(, flux of 93cm~'^s~^ or an event rate of 1.6 events (kton-year)~^ 
above the energy threshold (£77^ > 8.3 MeV). The present data also allows us to set 
more stringent limits on the diffuse supernova neutrino flux and on the annihilation 
rates for light dark matter particles. 

Subject headings: dark matter, ISM: supernova remnants, neutrinos, Sun: particle emis- 
sion 



1. Introduction 



Ambient electron a ntineutrinos {ves) of terrestrial origin include geoneutrinos ()Araki et al 



2OO5I : iBellini et al]l201o|) which have energ ies below ~3.4MeV, and man-made reactor antineu 



trinos (jCowan et al.l 1 19561 : IVogel et al 



1981 



y), which have energies below ~8MeV. Naturally pro- 
duced z/g's with higher energies must be of cosmic origin. The region above a few tens of MeV 
is dominated by atmospheric neutrinos generated from the decays of muons and pions produced 
by cosmic-ray interactions. There is a gap in energy between terrestrial and atmospheric neutri- 
nos where other "exotic" mechanisms to generate neutrinos could dominate. For antineutrinos 
in the energy region betwe en 8 MeV and 15 MeV, only diffuse neutrino flux from distant super- 
novae (iTotani &: Satdll995l) and exo tic generation mechanisms, e.g., co nversion of solar neutrinos 



into antineutrinos (lOkun et al.lll986l ) or light dark matter annihilation (jPalomares-Ruiz &: Pascoli 
20081 ) are thought to be possible. 



1.1. Solar Antineutrinos 



The nuclear fusion reaction produces most of the Sun's energy, and a portion of its energy 
is taken away by electron neutrinos. The data for all solar neutrino experiments are consistent 
with the prediction based on the Standard Solar Model including the flavor transition due to the 
neutrino oscillation. On the other hand, antineutrinos are produced by the f3~ decay of the natural 
radioactivity or the photo-fission of heavy isotopes in the Sun, but bo th flux contributions at Earth's 
surface are negligibly small relative to the terrestrial antineutrinos (IMalaney et al.lll990l ). 



However, if the neutrino has a non-zero magnetic moment it could be converted into an an- 
tineutrino in the strong solar m agnetic field. Thi s mechanism was originally p roposed as a solution 
to th e solar neutrino problem (|Okun et al.lll986l ) and was later revisited in (jAkhmedov &: Pulido 
2OO3I ). A two-step process takes place, the first step occurs deep inside the solar interior, where a 
fe converts into a via spin flavor precession. The 17^ then oscillates into an i7g while propagating 
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from the Sun to the Earth. The combined probabiHty for the two processes is 



P{UeL ^ l^eR) ^ 1-8 X IQ-^" sin^ 2^12 



■Q) 



lO-^^HB 10 kG 



(1) 



where Bt is the transverse solar magnetic field in the region of neutrino production, Rq is the solar 
radius, and ^ is the neutrino magnetic moment in Bohr magneton (/xb)- Very little is known about 
the magnitude of magnetic fields in the solar interior, but values up to 3 x lOj G are permitted 
based on Solar and Heliospheric Observatory observations (jCouvidat et al.ll2003l ). One can search 
for conversion of solar neutrinos because they have energies higher than terrestrial neutrinos. 
The present best limit for the probability of so lar neutrino-to-ant ineutrino conversion, from the 
Borexino experiment, is less than 1 .3 x 10~^ (Bellini et al.M201ll ) assuming an unoscillated ^B 
neutrino flux of 5.88 x 10®cm~^s~^ ( Serenelli 201ol ). The neutrino decay, predicting that a heavy 
neutrino mass eigenstate decays to a lighter an tineutrino mass eigenstate, is another possibility of 
antineutrino production (jBeacom &: Belli |2002| ). 



1.2. Diffuse Supernova Neutrino Flux 



A diffuse supernova neutrino background exists from distant core-collapse supernovae. The 
numerical calculation of supernova explosions predicts neutrino fluxes which are comparable to 
each other, and those energy spectra extend up to ~80 MeV. In typical hydrogen-rich detectors, 
the supernova neutrino signal is dominated by the z/g reaction due to the large cross section. Only 
upper limits for the diffuse supernova neutrino flux (DSNF) have been set by Super-Kamiok ande 
and the S udbury Neutrino Obse rvatory (SNO) for energies above 19.3 MeV (jMalek et al.ll2003l ) and 
22.9 MeV (jAharmim et al.ll2006l ) for I/g and respectively. Because the DSNF energy is lowered by 
redshift, the spectral shape is closely connected to the history of star f ormation. Various supernova 
and cosmological models predict different shapes (jAndo &: Satdl2004l ). while recent measurements 



of the core-co llapse supernova rate reduce the cosmological uncertainties (jHoriuchi et al 
Beacoinll2010l ). The model predictions should be confronted with data. 



1.3. Dark Matter Annihilation 



Dark matter self-annihilation into standard model particles can be searched experimentally. 
In light dark matter annihilation, most of the final state particles are kinematically forbidden, and 
only 77, e'^'e", and Vv are poss ible. However, th ere are strong limits on the branching ratio to 
the visible states, 77 and e~^e~ (jMack et al.ll2008l ). so the Vi' search bounds the total annihilation 
cross section. Because even a few MeV mass cold dark matter particles are non-relativistic, their 
annihilation should produce monoenergetic fluxes of neutrinos. Annihilation in the Milky Way halo 
should be the dominant source of such neutrinos, so redshift can be safely neglected. The Super- 
Kamiokande data provide the best limit on the rate of such annihilation for dark matter mass (m^) 
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above 15 MeV. Assu ming an annihilation model in the Galactic halo and equal annihilation into 
all neutrino flavors (|Palomares-Ruiz &: Pascolil 120081 ) . the limit can be translated into a velocity- 
dependent averaged cross section ((cr^v)) for dark matter particles. 



2. The KamLAND Experiment 

The KamLAND detector is located ~lkm under the peak of Mt. Ikenoyama (36.42°N, 
137.31°E) near Kamioka, Japan. The 2700 meters water equivalent (mwe) of vertical rock overbur- 
den reduces the cosmic-ray muon flux by almost flve orders of magnitude. A schematic diagram of 
KamLAND is shown in Figure [TJ The primary target volume consists of 1 kton of ultra-pure liquid 
scintillator (LS) contained in a 13 m diameter spherical balloon made of 135 /im thick transparent 
nylon EVOH (ethylene vinyl alcohol copolymer) composite film. The LS consists of 80% dode- 
cane and 20% pseudocumene (1,2,4-trimethylbenzene) by volume, and 1.36 it 0.03 g 1~^ of the fluor 
PPO (2,5-diphenyloxazole). A buffer comprising 57% isoparaffin and 43% dodecane oils by volume, 
which fills the region between the balloon and the surrounding 18 m diameter spherical stainless- 
steel outer vessel, shields the LS from external radiation. The specific gravity of the buffer oil is 
adjusted to be 0.04% lower than that of the LS. An array of photomultiplier tubes (PMTs) — 1325 
specially developed fast PMTs masked to 17 inch diaineter a nd 554 older 20 inch diameter PMTs 



reused from the Kamiokande experiment (jKume et al.lll983l ) — are mounted on the inner surface 



of the stainless steel, providing 34% photocathode coverage. This inner detector is shielded by a 
3.2 kton water-Cherenkov outer detector (OD). 

Electron antineutrinos are detected in KamLAND via the inverse beta-decay reaction, 

I7e + p — )• e+ + n. (2) 

This process has a delayed-coincidence (DC) event-pair signature which offers powerful background 
suppression. The energy deposited by the positron, which generates the DC pair's prompt event, is 
approximately related to the incident I7g energy by E^;^ ~ E'p + En + 0.8 MeV, where -E'p(= Tg+ + 
2772-6 ) is the sum of the kinetic energy and annihilation 7 energies, and En is the average 
neutron recoil energy which is below 1 MeV for the analysis Ep range, 7.5-30.0 MeV. Most of 
the neutron recoil energy is transferred to recoil protons, resulting in a negligible contribution of 
scintillation light due to the quenching effect. The delayed event in the DC pair is generated by a 
2.2 MeV 7-ray p roduced when th e neutron captures on a proton. The mean neutron capture time is 



(207.5 ± 2.8) yus (lAbe et al.ll2010l^ . The inv erse beta decay cross section is well approximated at the 



first order in 1/M (IVogel Beacomlll999l ). where M is the nucleon mass. The angular distribution 
of the positron emission is nearly isotropic, and unlike water Cherenkov detector, the scintillation 
light is also isotropic. As a result, the positron signal does not provide the incoming antineutrino 
source direction. Due to the extremely low cross section of antineutrinos, the Earth does not make 
a shadow extraterrestrial antineutrinos, and the detector has isotropic sensitivity. 



The detector is periodically calibrated with 7 sources deployed from a glove box installed at 
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the top of the chimney region. The radioactive sources are 60Co, 68Ge, 203Hg, 65zn, 24iAm9Be, 
^^^Cs, and ^^''Po^^C, providing energy cahbration up to ~8MeV along the central axis of the 
detector. In addition, capture of neutrons on hydrogen and carbon provides energy calibration 
throughout the entire sensitive volume. The visible energy in the detector is measured from the 
number of detected photoelectrons and is corrected for event position, detector non-uniformity, 
and scintillator nonlinearity from quenching and Cherenkov light production. The overall vertex 
reconstruction resolution is ~12cm / \J i?(MeV) and energy resolution is 6.4%/y^^(MeV). Energy 
reconstruction of positrons with > 7.5 MeV (i.e., Eqj^ > 8.3 MeV) is ve rified using tagg ed ^^B 
/3~-decays (r = 29.1ms and Q = 13.4 MeV) generated via muon spallation (|Abe et al.l 120101 ). 

From the studies of Bi-Po sequential decays, the effective equilibrium concentration of ^ss^j 
and ^^^Th in the LS are (2.2 ±0.3) x lO^^^g g"^ and (4.8 ±0.3) x lO'^'^g g~\ prior to the start of 
the LS purification campaign in 2007. From the singles energy spectrum fit, the concentration of 
^''K is (2.2 ± 0.2) X 10~^^ g g~^. Those radioactive impurities are negligible in this study relative 
to other backgrounds, such as muon spallation products and external 7-rays. The LS purification 
further reduced the radioactive impurities. 



3. Event Selection 

The present analysis includes data accumulated between 2002 March 5 and 2010 July 23, cor- 
responding to 2343 live-days. For the present search the following criteria were used: the prompt en- 
ergy is required to be 7.5 MeV < Ep < 30.0 MeV, and the delayed energy to be 1.8 MeV < E^ < 2.6 MeV; 
a fiducial volume cut of -R < 6m on both prompt and delayed events, a time correlation cut of 
0.5 fj,s < AT < 1000 fis, and a spatial correlation cut of AR < 1.6 m, this cut is driven by the mean 
free path of the capture-7 in the LS rather than the diffusion distance of the neutron. The in- 
efficiency caused by neutrons escaping from the LS is negligibly small. We limited the prompt 
energy window up to 30 MeV, because in the higher energy, there are background events from oil 
Cherenkov muons which are untagged by the OD due to a small inefficiency. Spallation cuts were 
used to reduce backgrounds from long-lived isotopes, such as ^Li (r = 257ms and Q = 13.6 MeV), 
that are generated by cosmic muons passing through the scintillator: a 2 ms veto is applied to the 
entire detector volume after a non-showering muon for both prompt and delayed events, a 2s veto 
is applied after a showering muon (i.e., muons depositing more than 3GeV of energy above their 
minimum ionizing contribution) or non-reconstructed mu on, while a 2 s 3 m radius cylindrical cut is 



applied around well-reconstructed non-showering muons (lAbe et al.ll2010l ) for delayed events. The 
overall selection efficiency of the candidates is 92%, which is evaluated from a Monte Carlo (MC) 
simulation. 
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4. Background Calculations 

4.1. Random Coincidences 

Two uncorrelated events in the detector may accidentally coincide in time, space, and energy so 
as to pass the Vf. selection cuts. To estimate the background contribution from random coincidences, 
events were selected with the appropriate prompt and delayed energies but in an out-of-time interval 
of 0.2 s to 1.2 s after the prompt event. This out-of-time window is 10^ times longer than the 
time interval used for the selection, providing a high-statistics background measurement. The 
time distribution between prompt and delayed events in the range between 0.2 s and 1.2 s shows no 
correlation between these events. The random coincidence background for the analysis is determined 
to be 0.22 ± 0.01 DC-pairs. 



4.2. Reactor Antineutrinos 



The location of the KamLAND detector was selected for the copious Ve, fl ux from 56 Japanese 

The reactor Vp flux 



2011 



nuclear power plants in order to study neutrino oscillation (jGando et al.l 
at KamLAND dominates all other Vf, sources for < 7.5 MeV. However, the tail of the reactor 
neutrino energy distribution extends to higher energies. The flux comes primarily from the beta 
decay of neutron-rich fragments produced in the fission of four isotopes: ^'^^U, ^'^^U, ^^^Pu, and 
^^^Pu. For each reactor the appropriate operational records including thermal power generation, 
fuel burn-up, shutdowns and fuel reload schedule were u sed to calculate the fission rates . The 



resul t ing i/f, spectrum was calculated using the model of (ISchreckenbach et al.lll985l : iHahn et al 



19891 : IVogel et al.lll98ll ) taking neutrino oscillation into account. The same methodology was used 
for previous reactor z7e analyses ar id showed excellent agreement over a wide energy range between 
expected and detected events (jGando et al.l l201ll ) . The total number of reactor candidates 
having Ep > 7.5 MeV is calculated to be 2.2 it 0.7 events, including a ~10% event-rate increase due 
to energy resolution. 



4.3. Radioactive Isotopes 



Cosmic-ray m uons interacting with carbon nuclei in the scintillator produce a variety of ra- 
dioactive isotopes (jAbe et al.ll2010l ). Two of these isotopes, '^He (r = 171.7ms and Q = 10.7 MeV) 
and ^Li (r = 257.2ms and Q = 13.6 MeV), have decay modes with electrons and neutrons in the 
final state. Such decays create DC pairs similar to inverse beta decay and therefore represent a 
background in the present study. 

The combination of a 2 s veto of the detector after showering muons and a 2 s 3 m radius 
cylindrical cut after non-showering muons significantly reduces the contribution of this background, 
but cannot eliminate it completely. The ^Li isotope, which has a higher end-point value, longer life 
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time and a higher production rate, generates the majority of these background events after cuts. 
To determine the contribution from this background we selected ^Li candidates using the same cuts 
that were used for the selection of the Vf, candidates, but the muon veto was not applied. The ^Li 
rate was evaluated from the distribution of the decay time relative to all previous muons, using a 
wider energy window of 0.9MeV< < 15.0 MeV to reduce statistical errors. For a 6m fiducial 
volume, 2074 it 49 events were found after showering muons and 454 it 31 events in a 3 m radius 
cylinder around the muon track after non-showering muons. Twenty percent of these events occur 
in the energy region of interest 7.5 MeV< Ep < 15.0 MeV. The 2 s cut reduces the ^Li background 
from showering muons to less than 0.2 events for Ep > 7.5 MeV. As non-showering muons occur with 
a relatively high frequency (0.2 Hz), to avoid the drastic loss of exposure which would accompany a 
2 s full-detector veto we take advantage of the fact that non-showering muons can be relatively well 
tracked in the LS and instead restrict the 2 s veto to a 3 m radius cylinder around the muon track. 
A 2 ms full- volume veto after all tagged muons is also used to suppress any spallation neutrons. To 
measure the efficiency of these cuts the distribution of neutron captures as a function of distance 
from the muon track was examined, and we found only 5.9% of neutrons survive the 3 m radius 
cylindrical cut. The resulting number of ^Li from spallation background with Ep > 7.5 MeV and 
surviving the 3 m radius cylindrical cut and time cut is 4.0 zt 0.3 events. The dead time introduced 
by all these cuts is 9.2%. 



4.4. Fast Neutrons 



Fast neutrons outside the inner detector may cause backgrounds in the fiducial volume. A fast 
neutron can scatter on protons or carbon nuclei in the LS producing a scintillation signal followed 
by a neutron capture signal, mimicking an coincidence. A MC simulation of fast neutrons reveals 
that the dominant background contribution is caused by muon-induced cosmogenic neutrons. A 
2 ms veto after OD-tagged muons mostly eliminates this background, while OD-untagged muons 
and the OD inefficiency cause a residual background. The MC-based study estimates 3.2 it 3.2 
fast neutrons remain in the data set, where a conservative unc ertainty of 100% for the simulated 
neutron production rate by muons is assumed (jAbe et al.ll2010l ). 



4.5. Atmospheric Neutrino Interactions 



Charged current (CC) and neutral current (NC) interactions of atmospheric neutrinos with 
carbon atoms in the KamLAND scintillator are the most significant source of background. At- 
mospheric neutrino spectra from (jHonda et al.l 120071 ). calculated specifically for the KamLAND 
location, were used to estimate the contribution from these backgrounds. 

The CC reactions by atmospheric V^s generate an irreducible background. The contribution 
from atmospheric Vg^s is estimated to be ^0.06 events in the energy window 7.5MeV< Ep < 
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30.0 MeV, which is dominated by the reaction on protons, bec ause the cross section fo r carbon nuclei 
is estimated to be at least one order of magnitude smaller (|Kim Cheounll2009l ). Atmospheric 
v^'s and z7^'s could react with both protons and carbon nuclei to produce muons and neutrons. 
The amount of detectable energy is shifted lower for such reactions because a large fraction of the 
neutrino's initial energy is expended to produce the muon. On the other hand, such reactions are 
followed by muon decay and therefore manifest themselves as a triple time correlation between 
the prompt event, muon decay and neutron capture. In the event selection, we found one triple 
coincidence event accompanied by a muon decay signal in a decay time interval of 0.5-10 fis. This 
event was excluded from the candida tes. To calculate the contribution from these reactions, the 
cross sections from (jAthar et al.ll2007l ) were employed. The resulting background levels for reactions 
with a neutron in the final state are listed in Table [TJ We estimate 4.0 it 0.9 events in total. 
The tagging efficiencies of the muon-decay-coincidence signature are calculated to be 78.6% and 
(77.5 lb 0.2)% from mean life times in carbon of positive and negative muons respectively. The 
residual background, including an untagged contribution from (7.1 it 1.4)% of negative muons 
which capture rather than decay, is 0.9 it 0.2 events. The observed rate of a muon decay signal 
is ~0.2 events (kton-year)~^ , which is c omparable to the b ackgrounds from "invisible" muons in 
Super-Kamiokande in this energy range (iMalek et al.ll2003l ). 

The most challenging background to estimate is that from the NC interactions of all neutrino 
species with carbon. In these reactions the neutrino transfers only a fraction of its energy to the 
final products. It can eject a neutron from the carbon nucleus, leaving it in an excited state with 
multiple decay modes. We used the following procedure to calculate the contribution from this 
background: we integrated the momen tum transfer from a neutrino to a quasi- free neutron over 



the e ntire atmospheric neutrino spectra (jHonda et al.ll2007l ) using cross sections from (lAhrens et al, 



19871 ). We then accounted for the neutron binding energies for P-shell (18.7 MeV) and S-shell 



(41.7 MeV) configurations and the corresponding shell populations. We also assumed that the 
neutron was remo ved from the carbon atorn , leaving it in an excited state. All de-excitation 
modes reported in (jKamyshkov &: Kolbell2003l ) were taken into account. For each final product we 
converted the particle energy to visible energy in the detector using an energy scale model that 
includes nonlinearities from scintillator quenching. Most of the outgoing neutrons have a kinetic 
energy less than 200 MeV, and the resulting visible energies are concentrated in the lower energy 
region, typically less than 100 MeV. The de-excitation of ^^C* is dominated by 2 MeV gamma-ray 
emission, which has little effect on the energy spectral shape. In the analysis energy window, the 
position separation of energy depositions between prompt and delayed signals is ~ 60 cm, which is 
larger than that for the thermal neutron case (~ 40 cm). We calculate the contribution from this 
background to be 16.4 events with an estimated systematic uncertainty of 29% which is driven by 
uncertainties in the atmospheric neutrino flux and the cross section of NC neutrino interactions; 
see Table [21 



We also attempted to estimate the NC background using t he NUANCE software tool (version 
3), which simulates neutrino interactions and related processes (ICasperil2002l ). However, we found 
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the code overestimates this background rate by a factor of ~2 relative to the above calculation in 
the energy region under study, mainly due to an inaccurate cross section for intra-nuclear nucleon 
re-scattering and an unexpected ~25 MeV positive offset of outgoing neutron energies. We therefore 
do not use the NUANCE-based estimation in this analysis. 



5. Data Interpretation 

We observe 25 events after the cuts described in Section [31 Figure [5] illustrates the distribution 
for the prompt event position, the delayed energy, the spatial correlation, and the time correlation. 
There are five two-neutron candidates which may be caused by ^^C(n, 2n)^^C reaction of fast 
neutrons. The estimated number of backgrounds for Vf. detection summarized in Table [3] is 26.9 it 
5.7 events in the prompt energy window 7.5 MeV < Ep < 30.0 MeV. Figure [3] shows the event 
distribution as a function of prompt energy. The data set presented here contains 16 times more 
statistics than the first KamLAND publication on this subject, allowing us to verify the expected 
background contribution in the analyzed energy window. The data are analyzed using an unbinned 
maximum likelihood fit to the event spectrum. The estimate for ^Li and reactor I7g are rather 
robust, on the other hand, reliable data for NC interactions in the energy range of interest do not 
exist and the method we used to calculate this background contribution has large uncertainties. To 
avoid possible bias from modeling in the NC background calculation, the normalization of the NC 
events is a free parameter in the spectral fits. 

From the unbinned maximum likelihood fit, the allowed region for the NC background and 
the conversion probability from i/p to Vf, is shown in Figure |4l assuming an unoscillated neu- 



trino flux of 5.94 X 10 cm s ( Pena-Garay Sz Serenellill2008l ). For the NC-floated normalization 



analysis, the upper limit for neutrino conversion is 5.3 x 10~^ at 90% CL, which corresponds to a 
solar I7g flux of 93 cm^^s^^ or an event rate of 1.6 events (kton-year)~^ above the energy thresh- 
old (£'77^ > 8.3 MeV; containing 29.5% o f the total neu trino flux). This limit is a factor 2.5 



improvement over the previous limit in (jBellini et al.ll201ll ) due to 24 times more exposure. For 



comparison, the rate analysis in the energy range 8.3 MeV < E-p^ < 15.0 MeV gives a slightly more 
stringent limit of 1.4 events (kton-year)~^, if we use all the constraints on the background estimates 
including the NC background (Table [3]) . The fitted NC background assuming zero solar Ve events 
is 14.8^*^5 4 events, which is in good agreement with the calculation (16.4 it 4.7 events). 

The probability for solar neutrino conversion can be predicted by the models of spin flavor 
precession and Mikheyev-Smirnov-Wolfenstein large mixing angle solution oscillations in the Sun. 
If the conversion model for ®B neutrinos of Equation ([1]) without the distortion of the ^B spectrum 
is assumed, we obtain the following limit on the product of the neutrino magnetic moment (/i) and 
the transverse solar magnetic field in the region of neutrino production (Bt)'- 



(3) 



- 11 - 



using the value of 34° for the mixing angle (|Gando et al.l 120111 ) . The current best limit on the 
neut rino magnetic m oment is from the GEMMA spectrometer, ^-p^ < 3.2 x 10^^^ fi^ at 90% 
CL (jBeda et alj |2010| ) . Lack of knowledge of the value of Bt limits KamLAND sensitivity to 
the neutrino magnetic moment. 

There data also test other potential V^, sources. Assuming an energy spectrum from the refer- 
ence model (lAndo &: Satdl2004ll. w hich is consistent with a recent model reducing the cosmological 
uncertainties (jHoriuchi et al.ll2009l ). we found an upper limit for the diffuse supernova Ve flux of 



139 cm^^s^^ at 90% CL in the analyzed energy range. This limit is weaker than our solar Uf, 
flux limit due to the strong anticorrelation between the signal and NC background events am- 
plified by the similarit y in their spectral shapes. This flux limit corresponds to about 36 times 
the model prediction (jAndo Satd 12004 ) . indicating poor statistical power in constraining the 
cosmological models using the current KamLAND data. The upper limit for the monochromatic 
z7e fl ux at each energy can be trans lated to a limit for the dark matter annihilation cross sec- 
tion (jPalomares-Ruiz &: Pascolil 120081 ). The dark matter annihilation limit varies weakly over the 
dark matter mass range due to limited statistics. We obtain {(Jav) < (1-3) xlO^^^cm^s"^ at 90% 
CL in the mass range 8.3 MeV < < 31.8 MeV, as shown in Figure O This is the most stringent 
constraint on the annihilation cross section below 15 MeV. 

Finally, we also present model-independent upper limits for Fg fluxes, as shown in Fig ured The 

limits are given at 90% CL based on the rate analysis using the Feldman-Cousins approach (jFeldman fc Cousins 



19981 ) with 1 MeV energy bins, including all the constraints on the background estimates in Table El 
The KamLAND data provide the best limits in the presented energy range 8.3 MeV < E^;^ < 18.3 MeV, 
owing to the efficient detection by the DC method and large exposure. Given that data are back- 
ground limited, mainly from the atmospheric neutrino NC interactions, accumulation of additional 
statistics is unlikely to improve this limit significantly. 

In conclusion, we report the spectrum of high-energy I7g candidates found in the KamLAND 
data set accumulated over more than eight years of detector operation. The live time exposure 
corresponds to 4.53 kton-year. In the energy range from 8.3 MeV to 31.8 MeV, no excess of Ve 
events over the expected background consisting of mostly atmospheric neutrino NC interactions, 
cosmogenically induced radioactivity, and reactor neutrinos were detected. The data allow signif- 
icantly improved limits on solar conversion probability, and on DSNF and annihilation cross 
section of dark matter below 15 MeV. The present level of background indicates limitations for 
future studies of z/g's in this energy range using KamLAND. 

While a better detector location could eliminate ^Li background and suppress reactor neutrino 
background, atmospheric neutrino NC interactions will continue to present significant challenges 
for next-generation large LS detectors. In water Cherenkov detectors, like Super-Kamiokande, 
the contribution from NC backgrounds is expected to be small, because the recoil protons by 
knocked-out neutrons from ^^O should be below the Cherenkov energy threshold. If gadolinium is 
added, Super-Kamiokande will overcome the problem of large backgrounds from solar neutrinos, 
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spallation products, and invisible muon decays, by the DC technique (jBeacom fc: VaginsI |2004| ) . 
and gain th e ability to detect the diffuse supernova neutrino signals. A possible future LENA 
experiment (IWurm et alj 120071 ) . which will have a detector with about 50 kton of LS, also a. ims to 
measure DSNF, reducing NC backgrounds by the ^^C tagging methods (jWurm et al.ll201ll ). For 
success, the detector design needs to be optimized to maximize the efficiency of the NC background 
rejection. 
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MEXT, Japan; and under the U.S. Department of Energy (DOE) grants no. DE-FG03-00ER41138, 
DE-AC02-05CH11231, and DE-FG02-01ER41166, as well as other DOE grants to individual in- 
stitutions. The reactor data are provided by courtesy of the following electric associations in 
Japan: Hokkaido, Tohoku, Tokyo, Hokuriku, Chubu, Kansai, Chugoku, Shikoku, and Kyushu 
Electric Power Companies, Japan Atomic Power Company, and Japan Atomic Energy Agency. 
The Kamioka Mining and Smelting Company has provided service for activities in the mine. 



APPENDIX 



The model-independent upper limits for Vf, fluxes provided for each energy may be useful to 
give an estimate of upper limits for various V^. sources. For example, one can easily test one's own 
model with a certain energy spectrum by an appropriate data integration. Table U] lists the 1 MeV 
binned upper limits shown in Figure [6l The binned is defined as 

v2 = "1 (4) 

where Vi is the model expectation for each energy bin, Ui is the KamLAND upper limit at 90% CL, 
and -v/2.71 is the conversion factor of limits from 90% CL to Icj CL. This binned analysis give 
the upper limit of the solar Vf. flux (< 87 cm^^s^^), and approximately reproduces the limit which 
is based on the unbinned maximum likelihood method including all background and systematic 
uncertainties. On the other hand, a limit for the diffuse supernova flux based on Equation ^ 
will be optimistic, because Table U] data include the constraints on the NC background estimate. 
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Fig. 1. — Schematic diagram of the KamLAND detector. 
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AR (m) AT (ms) 

Fig. 2. — Event distribution after all cuts: (a) prompt event position, (b) delayed energy, (c) spatial 
correlation, and (d) time correlation. The data are compared to the expected Vf, signal generated 
by the Monte Carlo simulation (lines). Dashed lines indicate the selection criteria for Vf. candidates. 
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Fig. 3. — Prompt energy distribution of selected Vf, candidates in 1 MeV bins together with the best- 
fit backgrounds (filled areas) and 90% CL upper limit of solar I7g's (red dashed line), which includes 
the detector response. The background histograms are cumulative and the solar Vf, histogram sits 
on top of the background total. 
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Fig. 4. — Allowed region, with the best-fit point (3.0 x 10~^, 14.5), for the unconstrained neutral 
current (NC) background and the probability of solar neutrino conversion from the KamLAND 
data in the energy range of 8.3 MeV < Ey^ < 31.8 MeV. The confidence level (CL) is shown for two 
degrees of freedom. The gray shaded region indicates the it la prediction from the NC background 
calculation. 
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Fig. 5. — Upper limits at 90% CL on the dark ma tter annihilation cross section fro m KamLAND 
(solid line) and Super-Kamiokande (dashed line) (jPalomares-Ruiz k. Pascolil l2008l ) . The shaded 
curve shows the natural scale of the annihilation cross section. 
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Fig. 6. — Model independent upper limits at 9 0% CL on the flux from KamLAND (blue line), 
Borexino (IBe lini et al.ll201ll) fr ed line), SNO (jAharmim et alJ l2004l ) (magenta line), and Super- 



Kamiokande (jGando et al.ll2003l ) (green line) . The 17^ flux limit depends on the event rate limit and 
the cross section o f Vp-P- The shaded curve shows the diffuse supernova flux for the reference 
model prediction (jAndo &: Satd |2004| ). Above the shown energy range, the Super-Kamiokande 

jerno va flux, owing to less 



data give more stringent upper limit of 1.2 cm ^s ^ for the diffuse su 



amount of backgrounds from muon spallation products (jMalek et al. 



20031 ). 
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Table 1. Calculated Backgrounds for Atmospheric Neutrino CC Interactions with Prompt 

Energy between 7.5 MeV and 30.0 MeV 



Reaction 


Number of Events 


Number of Untagged Events 


T'f, + p ^ + n 


2.1 


0.5 




0.7 


0.2 




0.4 


0.1 




0.4 


0.08 


T^ix+^^C^ /x+ + n+^Li+a 


0.4 


0.08 


V^+^^C^ H+ + 2n+iOB 


0.02 


0.005 


Total 


4.0 ± 0.9 


0.9 ±0.2 



Note. The numbers in the third column include the inefficiencies of the muon decay. 



Table 2. Calculated Backgrounds for Atmospheric Neutrino NC Interactions with Prompt 

Energy between 7.5 MeV and 30.0 MeV 



Reaction 


Number of Events 


z/(l7)+i2c^ + n+iiC+7 


13.2 




1.4 


u{u) + ^^C^ v{V) + n+^'Li+a + p 


1.4 


i/(i7)+12q_^ _^ n+'^Be+2j3 


0.3 


v{V)+^'^C^ u{v) + 2n+i°C 


0.1 


Total 


16.4 ±4.7 



Note. There is no muon in the final state, so the muon decay tagging is not useful unlike CC 

interactions. 
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Table 3. Summary of the Estimated Backgrounds with Prompt Energy between 7.5 MeV and 

30.0 MeV 



Background 


Number of Events 


Random coincidences 


0.22 ±0.01 


Reactor z/g 


2.2 ±0.7 


9Li 


4.0 ±0.3 


Fast neutron 


3.2 ±3.2 


Atmospheric v (CC) 


0.9 ±0.2 


Atmospheric v (NC) 


16.4 ±4.7 


Total 


26.9 ±5.7 



Table 4. Model Independent Upper Limit on the Ve Flux for Each Energy Bin from KamLAND, 

as Shown in Figure [6] 



Energy Range (MeV) 


Upper Limit at 90% CL (cm-^g-i) 


8.3-9.3 


56.2 


9.3-10.3 


67.1 


10.3-11.3 


63.8 


11.3-12.3 


14.0 


12.3-13.3 


25.8 


13.3-14.3 


32.2 


14.3-15.3 


18.9 


15.3-16.3 


24.1 


16.3-17.3 


14.5 


17.3-18.3 


13.0 



